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The problem of detecting specific features of microscopic dynamics in the macroscopic behavior 
of a many-degrees-of-freedom system is investigated by analyzing the position and momentum time 
series of a heavy impurity embedded in a chain of nearest-neighbor anharmonic Fermi-Pasta-Ulam 
oscillators. Results obtained in a previous work [M. Romero-Bastida, Phys. Rev. E 69, 056204 
(2004)] suggest that the impurity does not contribute significantly to the dynamics of the chain and 
can be considered as a probe for the dynamics of the system to which the impurity is coupled. The 
(r, t) entropy, which measures the amount of information generated by unit time at different scales 
r of time and r of the observable, is numerically computed by methods of nonlinear time-series 
analysis using the position and momentum signals of the heavy impurity for various values of the 
energy density e (energy per degree of freedom) of the system and some values of the impurity mass 
M. Results obtained from these two time series are compared and discussed. 
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I. INTRODUCTION 

The physically relevant question of how the underlying 
microscopic dynamics of many-particle systems is related 
to the observed macroscopic behavior is still not satisfac- 
torily answered. Since the main features of the latter 
are accounted for by the microscopic dynamics, which is 
chaotic for most systems of interest, it is natural to as- 
sume that some characteristics of the chaotic microscopic 
dynamics should be detected, in principle, at the macro- 
scopic level of description. Only then it would be possi- 
ble to assess the relevance of the microscopic dynamics 
on macroscopic behavior. 

Empirical evidence suggesting the existence of micro- 
scopic chaos on a molecular scale has been presented in 
an experiment on the position of a Brownian particle 
(BP) in a fluid The measurements were made at 
regular time intervals and the experimental time series 
data was then interpreted using standard techniques of 
time series analysis, suggesting a positive lower bound on 
the Kolmogorov-Sinai entropy per unit time h K3 , hence, 
microscopic chaos. However, a similar bound has been 
obtained with computer experiments on the nonchaotic 
Eherenfest wind-tree model where a single particle dif- 
fuses in a plane due to collisions with randomly placed, 
fixed, oriented square scatterers |2| • Later on it was found 
that the chaotic Lorentz model, which has circular scat- 
terers, has a diffusive behavior that is indistinguishable 
from the one exhibited by the Eherenfest model 3] . Fur- 
thermore, a class of one-dimensional maps has also been 
reported which present normal diffusivelike behavior in 
the absence of chaos, just as in the Eherenfest model 0. 
These results have rendered doubts about the conclusion 
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that microscopic chaos has been experimentally detected. 
The positive bound on h KS of Ref. pj has been attributed 
to the fmiteness in both spatial and temporal resolution 
and the limited amount of data points 0, ■ 

The results in all previous works have been obtained 
with only one of the BP dynamical variables, namely 
its position. However, there is a time scale in Brown- 
ian motion, not mentioned in the previous works, that 
is characterized by the relaxation time of the momen- 
tum autocorrelation function (MACF) of the BP. This 
function can be computed from the BP momentum time 
series. Significant variations in the BP momentum ocurr 
in this time scale, which is much smaller than the corre- 
sponding position time scale. Thus the BP momentum 
could, in principle, be used to probe the microscopic dy- 
namics of the system. But the type of the models used 
in Refs. 0- El lJ LJ precludes a study in the momentum 
space from the very beginning, and therefore a different 
system is needed to address the problem at hand. 

The Fermi-Pasta-Ulam (FPU) model, which is a one- 
dimensional chain of nearest-neighbor anharmonic oscil- 
lators, is a system that has been extensively studied over 
the past decades to test the way in which microscopic 
dynamics determines the macroscopic behavior. Now, if 
the FPU chain is coupled to a heavy impurity, the effects 
of the microscopic dynamics on the statistical behavior 
of the latter can be readily studied for a number of rea- 
sons. First, the position and momenta of all the oscil- 
lators and the heavy impurity are well defined variables. 
At variance with the Eherenfest and Lorentz models of 
Refs. in which the tracer particle interacts with the 

fixed scatterers through a discontinuous repulsive poten- 
tial, the system composed of the FPU chain coupled to 
the heavy impurity is defined in terms of a more realistic 
continuous attractive potential. Furthermore, this is a 
many-particle system in which the statistical behavior of 
the heavy impurity is produced by the interaction with 
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the oscillators of the chain, in sharp contrast with the 
models of Refs. 0, ||, which arc one-dimensional maps 
with only one dynamical variable, i.e. the position, de- 
fined in the model [a. 



II. THE MODEL AND ITS NUMERICAL 
INVESTIGATION 

In terms of dimensionless variables, the Hamiltonian 
of the model considered in this work is 



It has been known for some time that, for the FPU 
chain, a single parameter, namely its energy per degree 
of freedom e = E/N, controls its phase-space dynamics; 
that is, the system is weakly chaotic and thus has a dy- 
namical behavior that, in the timescales so far studied, is 
almost identical to that of a chain of harmonic oscillators 
for low e values, whereas it is strongly chaotic for high e 
values 7]. This behavior, for both low an high e values, 
is unaltered by the inclusion of a heavy impurity . The 
transition between these two dynamical regimes was first 
detected by means of a change in the scaling behavior of 
the largest Lyapunov exponent (LLE), which measures 
the exponential rate of divergence of two originally close 
trajectories in phase space, as e is varied from a low to a 
high value for the FPU model 7] and by a change in the 
scaling of the relaxation time of the MACF as a function 
of e for the FPU chain coupled to a heavy impurity @. 



In this work we will concentrate on the macroscopic 
manifestation of the transition from one dynamical 
regime to another, which can be readily studied in the 
macroscopic behavior of the heavy impurity embedded in 
a FPU chain. Our starting point comes from the observa- 
tion that the Kolmogorov-Sinai entropy h KS , which has 
been computed from the complete microscopic dynam- 
ics of this system as a function of e, displays the afore- 
mentioned transition in the scaling behavior of h KS 0. 
Since in the description of the microscopic dynamics the 
momentum of all the oscillators of the system is consid- 
ered and the dynamics of the chain is not affected by the 
presence of the heavy impurity, we can ask whether some 
evidence of the above mentioned transition between dy- 
namical regimes can be detected by applying the methods 
of nonlinear time scries analysis to the momentum time 
series of the heavy impurity. The present work is a first 
step to answer the above posed question. 
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with mi — 1 for i ^ and too = M; periodic boundary 
conditions are assumed (^(jv/2)+i = x -N/2)- The model, 
which will be named modified FPU (MFPU) model from 
now on, describes a system of one-dimensional TV coupled 
nonlinear oscillators of unit mass with nearest-neighbor 
interactions, displacements {x{\, momenta {pi} and a 
central oscillator (impurity) of mass M with displace- 
ment xq = X and momentum po = P. The value (3 = 0.1 
was used in the computation of all the numerical results 
hereafter reported. As initial conditions we choose the 
equilibrium value of the oscillators displacements, i.e. 
3^(0) = for i = -N/2, N/2. The momenta fe(0)} 
were drawn from a Maxwell-Boltzmann distribution at a 
temperature consistent with a given value of the energy 
density e, which was chosen in the range 0.01 < e < 100 
with M = 40, 60, 80, and 100. Then the 2(N + 1) equa- 
tions of motion were numerically integrated to obtain the 
time evolution of the system, whose state is represented 
by the variable T(t) = ({xi(t)}, {p l {t)}) G 3? 2<N+1) . After 
thermal equilibrium between the impurity and the FPU 
chain with N = 300 000 unit mass oscillators is attained, 
the behavior of the heavy impurity is almost identical to 
a Brownian motion for all e values studied 10], as can be 
inferred from the values of the diffusion coefficient and 
the exponential fit to the MACF @. So, we can consider 
the heavy impurity as a genuine BP. It is in this regime 
that the position {X(t a )} and momentum {P(t a )} time 
series of the heavy impurity, with t a = err {a = 1, . . . , n) 
and t = 1 being the natural time unit (discretization 
step), were recorded for all e and M values considered 
over a time interval of n = 2 x 10° natural time units. 



III. POSITION TIME SERIES ANALYSIS 



The plan of the paper is as follows. In Sec. II we 
briefly review the model and the methodology employed 
to obtain the position and momentum time series of the 
BP. In Sec. Ill we present the results for the position time 
series, which are consistent with those of Refs. Q,0|. In 
Sec. IV we apply the same methodology of the previous 
section to the momentum time series, with the result that 
there is a quantitative difference in the results obtained 
in the low and high e regimes, and hence evidence of the 
microscopic dynamics of the oscillator chain. In Sec. V 
we discuss the results so far obtained. Conclusions are 
given in Sec. VI. 



In Fig. n we show a segment of the complete position 
time series for e = 0.01 and 10 with M = 100, which 
correspond to the weakly chaotic (i.e. almost periodic) 
and one instance of the strongly chaotic regime, respec- 
tively. We can readily appreciate that both signals are 
practically indistinguishable (except for the scale, which 
is controlled by the e value), although they correspond to 
completely different dynamical regimes, as already men- 
tioned. Furthermore, some quasirecurrences can be ap- 
preciated in the complete data series (not shown). This 
fact is a consequence of the finite N value. Thus it can 
be inferred that the position of a BP is unsuitable as a 
probe to the microscopic dynamics of the chain. 
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FIG. 1: (a) Position time record for e — 0.01. The numerically 
computed diffusion coefficient is D m 0.0052. The length of 
the complete data set is 2 x 10 5 and the data were sampled 
with t = 1. (b) Position time record for e = 10 with a 
numerically computed diffusion constant of D ~ 33.65. Same 
total data set length and r value as in (a). In both cases 
M = 100. 

From the BP position time series alone the investiga- 
tion on the system is performed at a finite resolution r; 
moreover, the dimensionality of the phase space is lost, 
since X(t) 6 5ft. To cope with these limitations, the phase 
space is reconstructed by delay embedding technique 01 
in which, from the scalar time series {X(t a )}, a new vec- 
tor time series is defined as 

X^ n \t a ) = {X(t a ), X(t a + r), . . . , X(t a + (to - l)r)}, 

(2) 

being a portion of the discretized trajectory with step 
r and defined in 5ft m ; to is known as the embedding di- 
mension. Next, the space 5ft is partitioned using cells of 
length r, which defines a resolution for the X variable. 
The vector X( m )(i) is coded into a word of length to, 

X. (m \t a ) W m (r,t a ) = (i{r,t a ),...,i(r,t a +(m-l)T)), 

(3) 

where i(r, t a + jr) labels the cell in 5ft containing X(t a + 
jr). Under the hypothesis of stationarity the probabili- 
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FIG. 2: /i (m) (r,r) computed using 7 x 10 4 data points from 
the time series shown in Fig. The reported results are for 
an embedding dimension of m = 50. (a) e = 0.01 and (b) 
e = 10. The two straight lines show the corresponding D/r 2 
behavior. 

ties P(W m (r)) of the admissible words {W m (r)} are ob- 
tained from the time evolution of ~Xy n )(t a ). Then the 
(r, t) entropy per unit time h^ m \r,T), a generalization 
of h KS for variables measured with finite resolution |12| . 
is defined by 

fc("0(r,r) = -[H m+l {r,r) - H m (r,r)}, (4) 
r 

where H m is the TO-block entropy 

H m {r,r)=- J2 P(W m (r))lnP{W m (r)), (5) 

{W m {r)} 

which was computed by means of the Grassberger- 
Procaccia method 13j. The result for the e values re- 
ported in Fig. n is displayed in Fig. |3 

As can be readily seen, the shape of M m )(r,r) is very 
similar for the weakly chaotic case e = 0.01 and the 
strongly chaotic case e — 10. This result can be read- 
ily explained by the fact that the BP, in both dynamical 
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regimes, performs Brownian motion |8|- Now, a Brown- 
ian signal is a stationary Gaussian process characterized 
by a power spectrum S(u>) oc c<j~ 2 , which allows to ob- 
tain an explicit form for the (r, r) entropy in the limit 
r — > as h^ m \r, r) ~ D/r 2 , where D is the diffusion co- 
efficient This asymptotic behavior is also displayed 
in Fig. for each e value considered using the numeri- 
cally computed values of the diffusion coefficient in each 
dynamical regime Q ■ As in Refs. Q, 0, we consid- 
ered different delay times r to verify that the obtained 
result is indeed independent of the specific delay time 
value employed [&'). From the figure it is evident that 
this power law scaling for ft/ m )(r, r) is approximately ob- 
tained already for t = 1, with no noticeable differences as 
r increases. This result can be explained by the fact that 
the natural time unit r = 1 , employed as the lowest delay 
time, is very close to the value of the fastest period of the 
harmonic part of the chain, T m i n = 7r, which is the char- 
acteristic microscopic time scale of the system; thus the 
employed sampling time has an adequate resolution to 
probe the dynamics of the BP position. From the results 
displayed in Fig. |3 we can conclude that the stochastic 
behavior of the obtained signals is indeed independent of 
the underlying dynamics. In fact, as will be latter ex- 
plained, the mechanisms that explain the stochastic-like 
behavior for both e values are well known for both dy- 
namical regimes. Since our results are consistent with 
those of Refs. ^ ^ s cl ear that no microscopic infor- 

mation can be obtained from the position time series of 
the heavy impurity. 

IV. MOMENTUM TIME SERIES ANALYSIS 

Now we consider the momentum time series of the 
heavy impurity, which is displayed in Fig. for e = 0.01 
and 10. The first feature that stands out in compari- 
son with the corresponding position time series is that 
the time scale in which significant variations in the sig- 
nal value ocurr is much shorter than the corresponding 
time scale of the position variable. The second feature is 
that no quasirecurrences can be found by inspecting the 
complete time series (not shown), in sharp contrast with 
what happens with the complete position time series. Of 
course, for a long enough record quasirecurrences could 
be found, since N is finite, as already mentioned. How- 
ever, the fact that for the same data length the position 
variable shows quasirecurrences and the momentum time 
series does not clearly indicates that it may be possible 
to extract more dynamical information from this variable 
for both e values considered than from the corresponding 
time record of the BP position variable. 

As in the case of the position, a reconstruction of the 
phase space is performed from the BP momentum time 
series, which yields 

P (m) (*a) ={P(t a ),P(t a +T),...,P(t a + (m-l)T)}. (6) 

The (r, r) entropy per unit time for e = 0.01, obtained 
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FIG. 3: (a) BP momentum time record for e = 0.01. The 
length of the complete data set is 2 x 10 5 and the data were 
sampled with r = 1. (b) BP momentum time record for 
e = 10. Same total data set length and r value as in (a). 
In both cases M — 100 and the inset figures display the BP 
momentum time record over a larger time scale. 



from the vector time series P^(t a ), is displayed in Fig. 
01 for various m values. The most striking feature is that 
the (r, r) entropy has a well defined saturation value 
Zi'm^r, t) cons i ; which is reached asymptotically as 
the embedding dimension m increases, for a small reso- 
lution r. Furthermore, as can be readily seen, there is 
no need to take large delay time values in order to ob- 
tain the plateau value displayed. The finite value of the 
(r, t) entropy is consistent with the fact that the LLE, 
and thus the h KS entropy, are well defined quantities in 
the thermodynamic limit, both for the FPU chain 0] 
and for the MFPU model 9] . It is also clear that there is 
no way to obtain the aforementioned plateau value using 
the position time series alone, as shown in Fig. [3 The 
behavior of the (r, r) entropy for low r values, displayed 
in Fig. is a confirmation that the momentum of the 
BP is a better probe of the microscopic dynamics of the 
oscillator chain than its position. However, the satura- 
tion value hS m >{r,T) w 0.086 obtained for r < 1 is higher 
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FIG. 4: h (m) (r, r) computed using 7x 10 4 data points from the 
BP momentum time series shown in Fig. |3Ja). The reported 
results are for a delay time of r = 1 and various values of 
the embedding dimension m. The horizontal dashed line is 
the mean value obtained from the average of the curves with 
m > 40 for r values in the range 0.19 < r < 0.59. 

than the true microscopic value h KS k, 0.12 x 10 -4 0. 

In Fig. [SI we present the (r, r) entropy for different val- 
ues of the energy density e for a single delay time r = 1 
and m = 60 in all cases, since this is a representative 
value of the embedding dimension, as can be seen in Fig. 
0J We observe that the curves are separated into two dif- 
ferent and well-defined regions, depending on the e value. 
In all cases a plateau value can be denned when e < 1. 
On the contrary, when e > 1, it becomes increasingly 
difficult to find a region in which a plateau value can 
be observed. So, from this figure we can establish that 
e T ~ 1 is a threshold between two different regimes that 
are well characterized by the behavior of h^ m \r,T), as 
already explained. The observed difference in the behav- 
ior of the (r, r) entropy can be readily interpreted as a 
macroscopic manifestation of the structural changes that 
the phase-space undergoes when going from a low to a 
high e value. Further details will be given in the next 
section. 



V. DISCUSSION 

Our results shed some new light on the problem of 
detecting features of the underlying microscopic dynam- 
ics by means of macroscopic measurements. Since we 
have access to both the complete microscopic dynamics 
of the system and the results of data analysis of the BP 
position and momentum times series, we can make com- 
parisons between "exact" microscopic information and 
macroscopic "measurements" which are not possible to 
make if other models are employed, as already mentioned 
in the Introduction. For a deterministic chaotic system 
one has < h KS < oo (h KS = for regular motion), 
whereas for a random process h KS = oo. Now, consid- 
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FIG. 5: /i (ra) (r,r) computed using 5 x 10 4 data points from 
the BP momentum time series for different energy density 
values. The reported results are for an embedding dimension 
of m = 60, with M = 100 in all cases. 



ering the (r, r) entropy computed from the BP position 
time series, there is hardly a way of discerning from Fig. [2 
whether the data were originated by a an almost periodic 
or chaotic system, since the behavior of h^ m ' (r, r) is very 
similar in both cases. The problem can be traced to the 
embedding dimension to, since its value is well bellow 
the dimensionality 2(N + 1) = 600 002 of the underly- 
ing dynamical system. This is an intrinsic limitation of 
the nonlinear time series methods which has been pre- 
viously established 0, E| ■ The aforementioned problem 
persists, but is less severe, when the BP momentum time 
series is employed in the computation of h^ m \r, r), since 
a well defined plateau value is obtained. This result can 
be explained by the fact that the BP momentum varies 
on a smaller time scale than the corresponding BP posi- 
tion time scale, thus allowing a more detailed sampling 
of the phase space. The problem becomes the overes- 
timation of h KS when it is approximated by the (r, r) 
entropy computed both from the BP position and mo- 
mentum time series. The Kolmogorov-Sinai entropy is 
an upper bound to the h^ m '(r, r) value, and thus the re- 
lation h^ m \r, t) < h KS is not satisfied, as can be readily 
seen in Figs. |21 and 0] Since the length of both the BP 
position and momentum time series is the same we can 
conjecture that this overestimation is a consequence of 
the small value of the embedding dimension to. Finally, 
we observe that the rate of decrease in the plateau value 
of the (r, r) entropy as to increases becomes too small to 
justify any further increment in the embedding dimen- 
sion, which is computationally expensive. 

From Fig. it is clear that the saturation regime of 
the (r, t) entropy becomes increasingly harder to define 
as the e value is increased. Therefore two regions can be 
unambiguously defined by a transition from one type of 
behavior to another. From a macroscopic, i.e. statistical, 
perspective it is not at all clear what kind of mechanism 
could be held responsible of this change in the behav- 
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ior of the (r, r) entropy as e is increased. Furthermore, 
this transition sets in at a very precise threshold value 
e T « 1, which in Ref. [jj was called the strong stochas- 
ticity threshold (SST). This threshold indicates a tran- 
sition between two different regimes in the microscopic 
dynamics of the FPU chain: it is strongly chaotic and 
phase-space diffusion is fast when e > e T , whereas is 
is only weakly chaotic (i.e. almost periodic) and phase- 
space diffusion is slowed down when e < e T . Then, from a 
microscopic Hamiltonian perspective, the change of be- 
havior detected by means of the (r, r) entropy can be 
straightforwardly interpreted as a macroscopic manifes- 
tation of the aforementioned transition described by the 
SST. We can say that the difference in the behavior of 
the (r, t) entropy for low and high e values means that 
some distinction between different dynamical regimes has 
indeed been observed by nonlinear time-series analysis. 

Now, it is important to recall that, in the particular 
case of e = 0.01, the h KS entropy has a finite, albeit very 
small, value 0. The apparent dynamical randomness, 
i.e. finite (r, r) entropy value, of the Eherenfest model 
has been previously explained by the structural disorder 
of the model , which is the static analog of our ran- 
dom initial conditions. Thus it is quite plausible that the 
finite information amount carried in them is reflected in 
the finite (r, r) entropy value displayed in Fig. 0] which 
corresponds to the weakly chaotic regime e = 0.01. The 
increase in the (r, r) entropy of the MFPU model for high 
e values displayed in Fig. [3] can be explained by the fact 
that it has been known for some time that the dynamics 
of the FPU model, for high e values, mimics a random 
process. Thus the LLE scaling exponent, which charac- 
terizes its power-like dependence on e, can be obtained 
by means of a random matrix approximation 0. This 
dynamical randomness produces an increase in the (r, t) 
entropy relative to its plateau value for e = 0.01. The pic- 
ture that finally emerges is that, due to the small value 
of the embedding dimension m, a finite value of the (r, t) 
entropy will most likely be obtained when an observable 
of a high-dimensional system like the MFPU model is 
examined with nonlinear time scries analysis techniques, 
independently of the precise nature of its microscopic dy- 
namics, either regular or chaotic. Nevertheless, a distinc- 
tion can be made between those regimes by inspection of 
the low r behavior of the (r, r) entropy, which has a well 
defined plateau value for the weakly chaotic regime and 
a monotonically increasing value for the chaotic regime, 
ft will be a matter of future research to explore if this 
type of distinction can be made in other models. 

We would want to end this section with some final re- 
marks. The central argument of this paper is that the 
information obtained from a given system can be very 
different depending on the dynamical variable chosen to 
be studied by time-series analysis methods. Thus, we 
have found evidence of the microscopic dynamics of the 
FPU chain by employing the BP momentum instead of 
the position. Now, as pointed out in Ref. Q, the tran- 
sition from an almost periodic to a chaotic behavior in 



the FPU chain has been observed in a change in the scal- 
ing behavior of the relaxation time of the MACF of the 
BP coupled to the oscillator chain. This is a quantity 
that can be readily measured in an experimental setup 
by means of standard techniques, such as neutron scatter- 
ing [l^. Thus, it could be possible to explore the micro- 
scopic dynamics of more complicated systems by means 
of this variable, instead of using time-series analysis tech- 
niques. Within the scope of these latter techniques there 
remains the possibility to estimate the (r, r) entropy by 
observing faster processes such as Johnson thermal noise 
in an electrical circuit, for example. In this particular 
case the analog of the BP velocity would be the electri- 
cal current. In Brownian motion the stochastic behavior 
of the BP velocity is produced by the interaction of the 
BP with the fluid particles. In Johnson noise interac- 
tion between the conducting electrons and the thermally 
vibrating atomic lattice of the wire give rise to a tem- 
porally varying electromotive force that is the analog of 
the stochastic force that appears in the Langevin equa- 
tion of Brownian motion. If a suitable microscopic model 
of the random electromotive force could be developed, a 
theoretical analysis along the lines of this paper could be 
performed, with results readily verifiable experimentally, 
since the wire current is a directly measurable quantity. 
Of course, any other process defined by a Langevin-type 
equation in which a microscopic model of the correspond- 
ing fluctuating force could be provided would be adequate 
to explore the microscopic dynamics by means of macro- 
scopic measurements. 



VI. CONCLUSIONS 

In this work we have presented positive evidence, ob- 
tained from a systematic study employing time-series 
analysis methods, that certain qualitative information on 
the microscopic dynamics of a FPU chain can be detected 
by analyzing the momentum time-series of a heavy im- 
purity coupled to the chain. This detection was made 
possible by the fact that the BP momentum varies in a 
much shorter time scale than the BP position. Since the 
BP momentum obeys a Langevin equation, our results 
open the possibility that a similar methodology could be 
applied to other processes if a certain variable that obeys 
a Langevin-type equation could be properly defined and 
a suitable microscopic dynamical model to account for 
the macroscopic randomness could be provided. 
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